'is 


REPORT  DOCUMENTATION  PAGE 


la.  HE  ’ORT  SECURITY  CLASSIFICATION 


)ULE 


AD-A205  073 


8£R(S) 


6a.  NAME  OF  -’’ERFORMING  ORGANIZATION 
Virginia  Polytechnic  Insittute 
&  State  University  _  cs 


6b.  OFFICE  SYMBOL 
(If  applicable) 


6c  AOORESS  ry,  State,  and  ZIP  Code) 

Blacksburg,  VA  24061 


8a.  NAME  OF  FUNDING  /  SPONSORING 
ORGANIZATION 


■AfOSE- 


8b.  t  ICE  SYr>.  .  OL 

(l  , pplicable ) 


8'c.  AOORESS  (City,  State,  and  ZIP  Code) 
Building  410 

Bolling,  AFB  DC  203.  6448 


lb.  RESTRICTIVE  MARKINGS 


Form  .proved 
0 MB  N-  .0704-0188 


3,  DISTRIBUTION /AVAILABILITY  OF  REPORT 

Ai'0SR/Br  sossa-** 


m  m  m'i 


Boiling  A>'"3  DC 


5.  MONITORING  ORGANIZATION  REPC 

AFOSRTK. 


NUMBERS) 

8  9-  0  20 


7a.  NAME  O'  MONITORING  ORGANIZATION 
AFOSR 


7b.  ADDRESS  (City,  State,  and  ZIP  Code) 

Building  410 


Bolling,  AFB  DC  20332-644 


9.  PROCUREMENT  INSTRUMENT  IDENTIFICATI 
AFOSR  85-0250 


10.  SOURCE  OF  FUNOING  NUMBERS 


1 1 .  TITLE  C  -elude  Security  Classification) 
Algorithms  for  Nonlinear  equations 


PROGRAM 

PROJECT 

TASK 

WORK  L  . 

ELEMENT  NO. 

NO. 

NO. 

ACCESSION  .. 

6U02F 

2304 

A8 

- 1 

12.  PERSON*. 

iavne  T.  _ 

Tla.  TYPE  OF  REPORT 
FINAL 


THOR(S) 

.'atson 


13b.  TIME  COVEREO 
FROM  '  /  1  /  8  5 


14.  OATE  OF  REPORT  (Year,  Month,  Day) 

^8/30/Rf  5  October  1,  1988 


15.  PAGE  COUNT 
3 


IS.  SUPPLEMENTARY  NOTATI 


17. 


COSATl  C 


FIELD 

GROUP 

U8-GROUP 


18.  S  3JEC7  TERMS  ( Continue  on  reverse  if  necessary  and  identify  by  block  number) 


19.  A8STF  CT  (Continue  on  r  -rse  if  necessry  and  identify  by  block  number) 

’  ’.is  project  volved  research  in  three  areas:  mathematical  software,  globally  conver¬ 
gent  1 .motopy  meth  is,  and  1  brid  quasi-Newton  algorithms  for  la.rge  scale  structural  Cotimi- 
zatio;  The  homot  py  resea  ;h,  concerned  mainly  with  low  dimensional  ferociously  nonlinear - 
. probl  ms,  centered  on  prove .  ;  convergence  rierems,  devising  homotopy  curve  tracking  algoirthm 
and  d  velopment  of  the  mathematical  softv-  ;e  package  Homopack.  The  structural  optimization 
research  concerned  optimization  algorithms  'or  very  large- sparse  nonlinear  problems,  where 
maintaining  sparsity  is  absolutely  necessary  and  even  matrix  multiples  are  costly.  Structu’a 
optimization  and  equilibrium  configuration  computation  via  quasi-Newton  and  homotopy  tech¬ 
niques  require  an  entirely  different  techmnlogy  for  quasi-Newton  and  homotopy  algorithms,  us- 
lng  aalistic  test  problems  from  structural  mechanics.  ,h  ;  (let) 

0  0  r- 


20.  DISTRI81  N  /  AVAILABILITY  OF  A  STRACT 
Cj  UNCLASSi.TEO/UNLIMITEO  □  AME  AS  RPT. 
22a.  NAME  OF  RESPONSIBLE  INDIVID  AL 
__Dr^_£harles_Holland_  i _ 

DD  Form  1473,  JUN  86 


G  OTIC  USER! 


21.  ABSTRACT  SECURITY  CLASSIFICATION 

_  CAlW~9 _ 

22b.  TELEPHONE  (Include  Area  Code) 

(202)  767-  5025 

22c.  OFFICE  SYMBOL  1 

-JSL—  _ J 

“rev/oui  editions  art  obsolete,  * 


-  ANNUAL  REPORT  FOR  AFOSR  85-0250 
ALGORITHMS  FOR  NONLINEAR  EQUATIONS 

Layne  T.  Watson 
Department  of  Computer  Science 
Virginia  Polytechnic  Institute  &  State  University 
Blacksburg,  VA  24061 

October  1,  1988 


e-< 


This  project  involved  research  in  three  areas:  mathematical  software,  globally  convergent  ho- 
motopy  methods,  and  hybrid  quasi-Newton  algorithms  for  large  scale  structural  optimization.  The 
homotopy  research,  concerned  mainly  with  low  dimensional  ferociously  nonlinear  problems,  cen¬ 
tered  on  proving  convergence  theorems,  devising  homotopy  curve  tracking  algorithms,  and  devel¬ 
opment  of  the  mathematical  software  package  IIOMPACK.  The  structural  optimization  research 
concerned  optimization  algorithms  for  very  large  sparse  nonlinear  problems,  where  maintaining 
sparsity  is  absolutely  necessary  and  even  matrix  multiplies  are  costly.  Structural  optimization 
and  equilibrium  configuration  computation  via  quasi-Newton  and  homotopy  techniques  require  an 
entirely  different  technology  than  for  small  dense  problems.  The  goal  of  this  research  is  to  explore 
and  develop  sparse  matrix  technology  for  quasi-Newton  and  homotopy  algorithms,  using  realistic 
test  problems  from  structural  mechanics. 


Accomplishments  under  AFOSR. 85-0250. 

First  year:  9/1/85-8/30/86. 

1)  Steve  Billups’  MS  thesis  “An  augmented  Jacobian  matrix  algorithm  for  tracking  homotopy 
zero  curves”. 

2)  A  fluid  mechanics  application  “Slow  viscous  flow  in  a  syringe”. 

3)  A  SIAM  Review  survey  paper  “Numerical  linear  algebra  aspects  of  globally  convergent  homo¬ 
topy  methods”. 

4)  Final  testing  and  documentation  of  the  ODE  based  part  of  HOMPACK. 

5)  Final  testing  and  documentation  of  the  normal  flow  algorithm  in  HOMPACK. 

6)  Preliminary  work  on  an  augmented  Jacobian  matrix  homotopy  algorithm. 

7)  Preliminary  work  on  sparse  homotopy  and  hybrid  algorithms. 


Second  year:  9/1/86-8/30/87. 

1)  Completion  of  theory  for  an  augmented  Jacobian  matrix  homotopy  algorithm. 

2)  Final  testing  and  documentation  of  the  augmented  Jacobian  matrix  algorithm  in  HOMPACK. 

3)  Preliminary  work  on  theory  and  code  for  sparse  versions  of  all  the  algorithms  in  HOMPACK. 

4)  Proof  of  the  applicability  of  homotopy  methods  to  finite  element  approximations  of  nonlinear 
two- point  boundary  value  problems  (collaboration  with  Ridgway  Scott  from  the  University 
of  Michigan).  The  resulting  paper  “Solving  Galerkin  approximations  to  nonlinear  two-point 
boundary  value  problems  by  a  globally  convergent  homotopy  method”  will  appear  in  SIAM 
J.  Sci.  Stat.  Comput. 
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5)  Five  fluid  mechanics  papers- 

•  “Dissipative  waves  in  fluids  having  both  positive  and  negative  nonlinearity”, 

•  “Magnetohydrodynamic  flow  and  heat  transfer  about  a  rotating  disk  with  suction  and 
injection  at  the  disk  surface  ’, 

•  “Magnetohydrodynamic  flow  past  a  porous  rotating  disk  in  a  circular  magnetic  field”, 

•  “Magnetohydrodynamic  flow  between  a  solid  rotating  disk  and  a  porous  stationary  disk”, 

•  “Rotating  magnetohydrodynamic  flow  in  a  circular  magnetic  field”. 

6)  A  structural  optimization  paper:  “Two  point  constraint  approximation  in  structural  optimiza¬ 
tion”. 

7)  A  paper  on  complementarity  theory:  “A  geometric  problem  in  simplicial  cones  with  application 
to  linear  complementarity  problems”  (collaboration  with  K.  G.  Murty  from  the  University  of 
Michigan  and  L.  M.  Kelly  from  Michigan  State). 

S)  A  tutorial  survey  paper  “Globally  convergent  homotopy  algorithms:  a  tutorial”  to  appear  in 
Appi.  Math.  Comput. 

9)  Invited  to  lecture  on  homotopy  methods  at  the  1986  ODE  Conference  (Albuquerque,  NM), 
the  ASCE  Structures  Congress  (New  Orleans,  LA),  and  the  AMS  Regional  Meeting  (Logan, 
UT). 

10)  ACM  Trans.  Math.  Software  paper  and  algorithm  “HOMPACK:  A  suite  of  codes  for  globally 
convergent  homotopy  algorithms’7,  representing  the  culmination  of  over  10  years  of  research 
and  development.  HOMPACK  contains  over  13,000  lines  of  code,  is  part  of  Argonne  National 
Laboratory’s  netlib  software,  and  has  been  distributed  to  over  50  sites  in  the  U.S.,  Mexico, 
Canada,  England,  Germany,  Israel,  Norway,  Argentina,  Australia,  and  Ireland. 

Third  year:  9/1/87-9/30/88. 

1)  Final  testing  and  documentation  of  the  sparse  versions  of  the  algorithms  in  HOMPACK, 
which  are  based  on  a  new  preconditioned  conjugate  gradient  algorithm  developed  by  Kamat 
and  Watson. 

2)  Four  papers  on  parallel  homotopy  methods: 

•  “A  globally  convergent  parallel  algorithm  for  zeros  of  polynomial  systems”, 

•  “Message  length  effects  for  solving  polynomial  systems  on  a  hypercube”, 

•  “The  granularity  of  parallel  homotopy  algorithms  for  polynomial  systems  of  equations”, 

•  “Granularity  issues  for  solving  polynomial  systems  via  globally  convergent  algorithms  on 
a  hypercube”. 

3)  Two  papers  on  polynomial  systems,  reporting  significant  advances  in  the  algorithmic  efficiency 
of  solving  polynomial  systems  (collaboration  with  A.  Morgan  from  General  Motors  Research 
Laboratories  and  A.  Sommese  from  the  University  of  Notre  Dame): 

•  “The  mathematical  reduction  of  a  system  of  equations  for  a  heart  dipole  problem”, 

•  “Finding  all  isolated  solutions  to  polynomial  systems  using  HOMPACK”. 

4)  Two  mechanics  papers  (collaboration  with  C.  Y.  Wang  from  Michigan  State  University): 

•  “Spinning  of  a  liquid  film  from  an  accelerating  disc”, 
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•  “Free  rotation  of  a  circular  ring  with  an  unbalanced  mass”. 

5)  Two  applications  of  sparse  homotopy  methods  to  large  scale  structural  optimization: 

•  “Tracing  structural  optima  as  a  function  of  available  resources  by  a  homotopy  method”, 

•  “Design  of  laminated  plates  for  maximum  buckling  load”. 

6)  A  theoretical  analysis  of  various  homotopies  for  the  linear  complementarity  problem  (collab¬ 
oration  with  A.  Poore  from  Colorado  State  University).  The  resulting  paper  “Continuous 
homotopies  for  the  linear  complementarity  problem”  will  appear  in  SIAM  J.  Matrix  Anal. 
Appl. 

7)  A  theoretical  analysis  of  quasi-Newton  updates  for  underdetermined  systems  of  equations 
(collaboration  with  H.  Walker  from  Utah  State  University),  “Least  change  secant  update 
methods  for  underdetermined  systems”. 

8)  Invited  to  lecture  on  parallel  algorithms  at  the  NATO  Advanced  Study  Institute  on  Numerical 
Linear  Algebra,  Digital  Signal  Processing,  and  Parallel  Algorithms  (Leuven,  Belgium),  and  on 
homotopy  methods  at  the  ORSA  CSTS  conference  (Williamsburg,  VA)  and  the  ORSA/TIMS 
Joint  National  Meeting  (Vancouver,  British  Columbia). 
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